Numerical solution of KKT systems in PDE-constrained optimization problems via the affine scaling trust-region approach

A recently proposed trust-region approach for bound-constrained nonlinear equations is applied to the KKT systems arising from the discretization of a class of PDE-constrained optimization problems. Two different implementations are developed that take into account the large dimension and the special structure of the problems. The linear algebra phase is analyzed considering the possibility of solving the arising linear systems by either direct methods or short-recurrence iterative linear solvers. Viability of the approach is proved through several numerical experiments on large KKT systems arising from the discretization of control problems.